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Abstract 

The form factors for the transitions ir"f, 777, 7/7 and t/ c 7 are ana- 
lyzed within the modified perturbative approach in which quark trans- 
verse degrees of freedom are retained. The results for the form factors 
are compared to experiment in detail. As compared to previous calcu- 
lations within the same approach only little modifications of the me- 
son distribution amplitudes are required in general in order to achieve 
agreement with experiment. Only for the irj form factor a strong 
contribution from the second Gegenbauer term is found. It also com- 
mented on the case of two virtual photons and on the transition form 
factors in the time-like region. 
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1 Introduction 



The simplest exclusive observable is the ir^y transition form factor, F 7r7 . It has 
been shown [1] that its behavior for large photon virtuality, Q 2 , is determined 
by the operator product expansion of the product of two electromagnetic 
currents near the light cone. The only soft physics information required in 
the calculation of the form factor is the pion wave function. It is however 
theoretically not understood what large Q 2 means or, in other words, at 



which value of Q 2 this collinear factorization approach can be applied. This 
is still to be decided by comparison with experiment. 

In 1995 the CLEO collaboration [2] showed the first data on the 7T7 tran- 
sition form factor at fairly large values of Q 2 , actually up to 8 GeV 2 . These 
data provoked an enormous theoretical activity. It seemed that collinear 
QCD provides a large contribution to the transition form factor. The data 
lie only about 20% below the asymptotic limit of that form factor, namely 
Q 2 F^y ->■ y/2f n where f n (= 131 MeV) is the pion's decay constant. Next-to- 
leading order (NLO) corrections [3] account for about a half of the difference 
between the asymptotic results and the data. For the remaining difference 
several explanations were proposed: a pion distribution amplitude slightly 
different from the asymptotic one ($as — Qx(l — x)), low renormalization 
scales which enhance the NLO corrections, power corrections or quark trans- 
verse momenta to mention a few mechanisms. Since all these mechanisms 
have to provide only little effects the ny transition form factor was believed 
to be the theoretically best understood exclusive observable. 

Recently, this believe has been ruined. The BaBar collaboration [4] has 
measured the transition form factors up to about 35 GeV 2 . While these 
data agree with the CLEO data [2, 5] below 8 GeV 2 they reveal an expect- 
edly large rise towards large Q 2 . This behavior is in dramatic conflict with 
dimensional scaling and turned previous calculations obsolete. As the CLEO 
data in 1995 the BaBar data renewed the interest in this quantity and many 
paper papers have already been devoted to its theoretical analysis, e.g. [6] - 
[13]. 

Here, in this article, it is proposed to employ the modified perturbative 
approach (MPA) invented by Sterman and collaborators [14, 15]. This ap- 
proach which bases on factorization, has been used before in the analysis 
of the CLEO data [16, 17]. It will be shown below that a good fit to the BaBar 
data can be achieved with this approach, only a few higher order terms in the 
Gegenbauer representation of the distribution amplitudes have to be taken 
into account now. Arguments will also be given why in the Q 2 range covered 
by the CLEO data the simple asymptotic distribution amplitude suffices for 
a fair fit to the data. In the next section the MPA will be described in some 
detail and its properties discussed. Next, in Sect. 3, the actual analysis of 
the data on the ivy form factor will be presented and the results compared 
with other theoretical approaches. Sect. 4 is devoted to the analysis of the 777 
and 7/7 transition form factors which have been measured by CLEO [5] and 
the L3 collaboration [18] previously and by BaBar [19] recently. The Babar 
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data [20] on the ?y c 7 transition form factor will be analyzed in Sect. 5 briefly. 
Comments on the case of two virtual photons and on the form factors in the 
time-like region will be presented in the respective Sects. 6 and 7 before the 
papers ends with a summary (Sect. 8). Details of the Sudakov factor will be 
presented in an Appendix. 

2 The modified perturbative approach 

An alternative to the usual collinear factorization is the transverse-momentum 
(kjj-factorization. For hard exclusive processes this type of factorization has 
been proposed by Sterman and collaborators [14, 15]. Further arguments for 
the validity of k^-factorization are given in [21] for the case of the transi- 
tion form factors. Nevertheless a rigorous proof of kj_-factorization does not 
exist as yet. The basic idea of kj_-factorization is to retain the quark trans- 
verse degrees of freedom in the hard scattering. This however implies that 
quarks and antiquarks are pulled apart in the transverse-configuration or 
impact-parameter space. The separation of color sources is accompanied by 
the radiation of gluons. Based on previous work by Collins et al [22, 23, 24] 
the corrections to the hard scattering process due to gluon radiation have 
been calculated in Ref. [14] in axial gauge using resummation techniques and 
having recourse to the renormalization group. These radiative corrections 
comprising re-summed leading and next-to-leading logarithms which are not 
taken into account by the usual QCD evolution, are presented in the form of 
a Sudakov factor in the impact-parameter plane. 

The k^-factorization combined with a Sudakov factor is termed the MPA. 
Its advantage is that the end-point regions where one of the parton momen- 
tum fractions tends to zero, are strongly damped. Large contributions to 
the form factor accumulated in the soft end-point regions would render the 
use of perturbation theory inconsistent. As has been pointed out by Isgur 
and Llewellyn- Smith [25] this is frequently the case in collinear factorization 
at experimentally accessible values of momentum transfer, typically a few 
GeV 2 . Another advantage of the MPA is that the renormalization scale, 
fj,R, can be chosen to be momentum-fraction dependent in order to elimi- 
nate large logarithms from higher-order perturbative corrections. Eventual 
a s -singularities are canceled by the Sudakov factor without introducing ad- 
ditional ad hoc cut-off parameters as for instance a parton mass. Thus, the 
MPA provides well-defined expressions for form factors and amplitudes of 
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Figure 1: Lowest order Feynman graphs for the P7 transition form factor. 

hard scattering processes: the perturbative contributions can be calculated 
in a self-consistent way, even for momentum transfers as low as a few GeV 2 . 

According to [16, 17] the form factor for a transition from a photon to a 
pseudoscalar meson (P) reads 

F Pl (Q 2 ) = J dx^ $ P (x, -b, n F )f§(x, b, Q, n^e-a^M****) , (1) 

within the MPA. Since the Sudakov exponent S is given in the impact- 
parameter (b) space, see App. A, it is convenient to work in that space which 
is canonically conjugated to the k^-space. In the convolution formula (1) T# 
is the Fourier transform of the momentum-space hard scattering amplitude 
Th evaluated to lowest order perturbative QCD from the Feynman graphs 
displayed in Fig. 1 but taking into account quark transverse momenta 



WQCf 
xQ 2 + k 2 A 



T£=^6r, (2) 



Here Cp is a charge factor. For instance, for the pion it reads C n = (e 2 — 
e|) / y/2 where e a denotes the charge of a flavor-a quark in units of the positron 
charge. The Fourier transform of (2) is 

fP = 2y/ECp (3) 

7T 

where the Fourier transform is defined by 

/(b) = j^y 2 J d 2 k ± exp [-i b ■ kj/(k ± ) . (4) 

The function Kq denotes the modified Bessel function of order zero. 

Another item in (1) is \l/p, the light-cone wave function of the meson P 
in the impact-parameter space. In the original version of the MPA [15] this 



wave function is assumed to be just the meson distribution amplitude $p. As 
argued in [15] the Sudakov factor e~ s can be viewed as the perturbatively 
generated transverse part of the wave function. For low and intermediate 
values of Q 2 , however, the non-perturbative intrinsic b- or k^-dependence of 
the light-cone wave function cannot be ignored as has been pointed out in 
[26] . The inclusion of the transverse size of the meson extends considerably 
the region in which the perturbative contribution to the form factor can be 
calculated. As in [17, 26] the wave function is parameterized in the form 



fp r xxh 

^ P (x,b,fi F ) = 2tt-j=§p(x,{i f ) exp 



Aa 2 P 



(5) 



The distribution amplitude, $p, possesses a Gegenbauer expansion 



<S> P (x,fi F ) = $ AS 



'^4X< n Cll\2x-l) 

n=2,4,.. « s (/! r 



1+ J2 °n(/ i o)(- 



(6) 



where the evolution of the expansion parameters a n from an initial scale, 
/i , to the factorization scale, /ip, is controlled by the anomalous dimensions 
(n > 0) 

7+2 = 7 + £L (2n + 5)(n* + 5n + 5) () 
7n+2 ' 7n fa (n + l)(n + 2)(n + 3)(n + 4)- [< ' 

Here, 70 = 0, C P = 4/3 and /3 = 11 — 2/3nj, rif is the number of active 
flavors. For the cases of interest in this work the distribution amplitude is 
symmetric under the replacement of the momentum fraction x by x = 1 — x. 
This symmetry is already taken into account in (2). The transverse size 
parameter op being related to the r.m.s transverse momentum by (with a = 

1) 



<k5 



67T 2 / 2 



TV 2 



■2 



dxx 2 x 2 E [al{CTf + 2a n a n+ ,Cl' 2 C^ 2 



5P, 



<n 



n=0,2,- 



1 - -a 2 (l - 2a 2 ) - — a 4 (a 2 - — a 4 ) + . . . 



(8) 



For later use the probability of the meson's valence Fock state is also quoted: 
d 2 b , s , „ o „ o -sr^ o (n + l)(n + 2) 



p«= /^^|^p(x,b)| 2 = 3tt 2 /I4 E « : 



n=0,2,- 



2n + 3 



(9) 
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The requirement P q g < 1 leads to bounds on the Gegenbauer coefficients or 
strictly speaking, on the products a n o P 



1 2n + 3 

a n a P <—-J— — — - — — -. 10 
7i f P y 3(n + l)(n + 2) 

For instance, if ap = 0.5 (1.0) GeV -1 at the initial scale /i one finds 

a 2 M < 2.14 (1.07) , a 4 M < 1.70 (0.85) . (11) 

The Sudakov exponent S in (1) which comprises the characteristic dou- 
ble logarithms produced by overlapping collinear and soft divergencies for 
massless quarks, is given in App. A. The impact parameter b which repre- 
sents the transverse separation of quark and antiquark, acts as an infrared 
(IR) cut-off parameter 1 [22, 23, 24]. Thus, 1/6 in the Sudakov exponent 
marks the interface between the non-perturbative soft momenta which are 
implicitly accounted for in the hadron wave function, and the contributions 
from soft gluons, incorporated in a perturbative way in the Sudakov factor. 
Obviously, the IR cut-off serves at the same time as the gliding factorization 
scale 

Hf = 1/6 (12) 

to be used in the evolution of the wave function. In accord with this interpre- 
tation the entire Sudakov factor is continued to zero whenever 6 > 1Aq CD . 
In this large-6 region the wavelength of the radiative gluon is larger than 
1/Aqcd- Because of the color neutrality of the hadron such gluons cannot 
resolve the hadron's quark distribution; hence radiation is damped. Soft glu- 
ons with wavelength larger than 1/Aqcd are therefore to be excluded from 
perturbation theory; they have to be absorbed into the soft wave function. 

Radiative corrections with wavelengths between the IR cut-off and the 
upper limit \/2/^Q yield to suppression through the Sudakov factor as is 
evident from the integration limits in (A. 2) (£ is either x or 1 — x). Gluons 
with even shorter wavelengths are regarded as hard ones which are considered 
as higher-order perturbative corrections of the hard scattering and, hence, 
are not part of the Sudakov factor. For that reason, the Sudakov function 



1 For a more complicated system like the proton's electromagnetic form factor, there are 
several b's. In order to cancel the a s singularities the b's have to be chosen appropriately 
[27]. 
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Figure 2: The exponential of the Sudakov function s(£i,bi,Q) vs £j = £ and 
6 = 6 for Q = 30A QCD . The hatched area indicates the hard scattering 
region. The plot is taken from [27]. 

s(£,&, Q) which is defined in (A. 6) and is related to S according to (A.l), is 
set equal to zero [15] whenever 

< 13 > 

In Fig. 2 the exponential of the Sudakov function exp[— s(£, b, Q)\ for Q = 
30 Aqcd is displayed. The properties of the Sudakov function lead to an 
asymptotic damping of any contribution except those from configurations 
with small quark- ant iquark separations, i.e. for InQ 2 — > oo the limiting be- 
havior of the transition from factors in collinear factorization emerges, for 
instance, F 7r7 — >■ y/2f T /Q 2 . 

In a NLO calculation of the P7 form factor or, say, in the case of the 
electromagnetic pion form factor a momentum-fraction dependent renormal- 
ization scale leads to a singular a s in the limit hr — > Aqcd- This singularity 
is canceled in the MPA: whenever a s tends to infinity the Sudakov factor e~ s 
rapidly decreases to zero. As one may suspect from Fig. 2 this does not seem 
to be the case in the region determined by 6Aq CD — > 1 and simultaneously 
£ < a/2Aqcd/Q, where exp[— s(£, b, Q)\ is fixed to unity. However, in this 
region the other Sudakov function in (A.l), namely s(l — £,b,Q), provides 
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the required suppression. The logarithmic singularities [In (1/(6Aq CD ) 2 ]~ 771 
which arise from the evolution of the wave function and which become worse 
with increasing Gegenbauer index are analogously cancelled by the Sudakov 
factor. 

In analogy to the case of the pion's electromagnetic form factor [15] the 
maximum of the longitudinal scale appearing in the scattering amplitude (2) 
and the transverse scale 

Hu = max(y/xQ, 1/6) , (14) 

is chosen as the renormalization scale [16, 17]. Although to lowest order there 
is no a s in the hard scattering amplitude for the P7 transition form factor, 
it nevertheless depends on /ir. Indeed, as discussed above, the Sudakov 
factor comprises the gluonic radiative corrections for scales between 1/6 and 
£Q/V2. Hence, the latter scale specifies the onset of the hard scattering 
regime. 

Inserting the Gegenbauer decomposition (6) of the meson distribution 
amplitude into (1) and integrating term by term one can write the P7 tran- 
sition form factor as 

Q 2 F P , = QCpf P ]T aMC n (Q 2 ,fi ,a P ) (15) 

n=0,2,4... 

The functions C n incorporate the integrals over the product of the n-th Gegen- 
bauer component of the meson wave function, the hard scattering amplitude 
and the Sudakov factor as well as the change of the Gegenbauer coefficients 
with the factorization scale. The evolution is worked out with the 1-loop 
expression for a s using Aqcd = 0.181 GeV and four flavors. This value 
of Aqcd is used throughout the paper except stated otherwise. In (15) /xq 
merely acts as the scale at which the Gegenbauer coefficients of the distri- 
bution amplitude are quoted; the transition form factor is independent of it. 
Numerical results for the ratios C n /C obtained under neglect of the intrinsic 
transverse momentum (i.e. for op — > 00), are shown in Fig. 3. A remarkable 
property of the MPA is to be observed from this figure: For n > and low 
Q 2 the Gegenbauer terms C n are suppressed as compared to the lowest one, 
Co. The strength of the suppression grows with the Gegenbauer index. A 
closer inspection of the functions C n in (15) reveals that the Sudakov fac- 
tor in conjunction with the hard scattering amplitude provides a series of 
power suppressed terms which come from the region of soft quark momenta 
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0) and grow with the Gegenbauer index [28]. With increasing 
Q 2 the higher Gegenbauer terms become gradually more important. At very 
large Q 2 however the evolution of the expansion coefficients again suppresses 
all C n except Co- Asymptotically, one has Co = 1 and C n = for n > and 
the asymptotic limit of the transition form factor emerges. The behavior of 
the C n in the MPA is very different from that in the collinear factorization 
approach. To LO for instance one has 

C coii = 1; C coii = [as{flF y as{fio)V n . (16) 




1 2 3 5 10 20 30 50 100 300 

Q 2 [GeV 2 ] 



Figure 3: The relative strength of the contributions from the n-th Gegenbauer 
term to the P7 form factor, scaled by the n = term, at the scale 1.0 GeV. 
Intrinsic is omitted. (Colors online) 

The intrinsic transverse momentum dependence of the wave function, 
on the other hand, also provides a series of power suppressed terms but 
these come from all x and do not grow with n [28]. These power corrections 
superimposes those generated by the Sudakov factor in combination with the 
hard scattering amplitude. For small values of the parameter ap the power 
corrections from the intrinsic transverse momentum dependence are rather 
strong and weaken the relative supppression of the higher Gegenbauer terms, 
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in particular for small n. In exchange for this the Cq term itself increases more 
sharply with rising Q 2 towards the asymptotic limit Co = 1. 

This property of the MPA explains why the CLEO data on the ir / ~f transi- 
tion form factor [5] are well described by the asymptotic distribution ampli- 
tude as shown in [17] (with a n = 0.861 GeV" 1 , A QC d = 0.20 MeV), see Fig. 
4. With the BaBar data [4] at disposal which extend to much larger values of 
Q 2 and do exceed the asymptotic limit y/2f w , higher Gegenbauer terms can 
no more be ignored; they are now required for a successful description of the 
transition form factors. What can be learned about the higher Gegenbauer 
terms from the BaBar data will be discussed in the next sections. 

3 Confronting with the BaBar data on tt7 
transitions 

Now, having specified all details of the MPA, one can analyze the 7T7 form 
factor by inserting (5) and (6) into (1) and fitting the Gegenbauer coefficients 
to the CLEO [5] and BaBar [4] data. From a detailed examination of the 
data it becomes apparent that besides the transverse size parameter only 
one Gegenbauer coefficient can safely be determined. If more coefficients are 
freed the fits become unstable. The coefficients acquire unphysically large 
absolute values between 1 and 10 (often in conflict with the bound (10)) and 
with alternating signs leading to strong compensations among the various 
terms. The reason for this fact is that the data show the tendency of a 
sharp increase of the slope in the vicinity of Q 2 — 10 GeV 2 while theory for 
a n > provides a slowly decreasing slope. With large absolute values and 
alternating signs of the coefficients a slightly increasing slope is produced for 
a finite range of Q 2 . 

Two fits are ultimately performed: for the first one the Gegenbauer coeffi- 
cient a 2 is fitted; for the second one a 4 is freed but a 2 is fixed at the face result 
of a recent lattice QCD calculation [29]. Evolved to the scale /i = 2 GeV 
with A QCD = 181 MeV, the lattice result [29] is a 2 (fi ) = 0.201 ± 0.114. All 
other Gegenbauer coefficients are assumed to be negligible in the fits. The 
results of the two fits are 2 

a* = 0.40 ±0.06 GeV" 1 , a 2 (^ ) = 0.22 ± 0.06 , x 2 = 34.1, (17) 

2 The Gegenbauer coefficient a n at a scale \i\ is obtained from the value quoted in (17) 
or (18) by multiplying it with the factor [In (MoMqcd)/ l n (^iMqcd)] 7 "- 
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Figure 4: The scaled 717 transition form factor versus Q 2 evaluated from fit 
(18) (solid line). The dashed line represents the result given in [17] which is 
obtained from the asymptotic distribution amplitude and a n = 0.861 GeV -1 
given. The dotted line is obtained from collinear factorization to NLO accu- 
racy. Data taken from [4, 5]. (Colors online) 

and with a 2 (/x ) = 0.20 [29], 

a n = 0.40 ±0.06 GeV -1 , a 4 (// ) = 0.01 ± 0.06, \ 2 = 34.2 , (18) 

The fits are stable and have a sharp \ 2 minimum. The minimal values of % 2 
for the two fits are reasonable given that 28 data points for Q 2 > 2.3 GeV 2 are 
included in the fits. The fit (17) yields a value for a-i that is in good agreement 
with the lattice result [29] within errors. The fit (18) is shown in Fig. 4; the 
fit (17) is practically indistinguishable from (18). At the largest values of Q 2 
the resultant form factor seems to be a bit small as compared to the BaBar 
data. Partially responsible for this fact are the large fluctuations the BaBar 
data exhibit; the fit compromises between all the data. At very large values 
of Q 2 , not shown in Fig. 4, the theoretical result for the form factor, scaled 
by Q 2 , flattens off with a broad maximum at about 200 GeV 2 and decreases 
subsequently towards its asymptotic value of v^/tt- The reason for the small 
value of ov in the present fits can easily be understood. The C2 and C4 terms 
demanded by the BaBar data [4] at large Q 2 , also contribute at low Q 2 , 
although to a lesser extent. Therefore, the Co-term with a w = 0.861 GeV -1 
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which fits nicely the low Q 2 data, see Fig. 4, is to be reduced. The only 
parameter available for this is however a n . 

In Ref. [30] it is argued that the decay 7r° — > 77 fixes the wave function 
at = integrated over x. For the wave function (5) this constraint can 
be turned into a result for the transverse size parameter [31] 

ra=0,2,4... 

For the asymptotic distribution amplitude this constraint leads to the value 
a w = 0.861 GeV -1 used in [17]. The parameters given in (17), (18) violate 
(19) plainly. Since it is however unclear at which scale the constraint (19) 
holds the significance of this violation cannot be judged. 

Also shown in Fig. 4 is a typical result of the collinear factorization ap- 
proach to NLO accuracy (taking /ip = lir = Q and, as an example, the 
Gegenbauer coefficients of the fit (18)). Apparently the shape of that result 
is in conflict with experiment. A change of the values of the Gegenbauer 
coefficients or the addition of further coefficients with the proviso that large 
negative coefficients are excluded, does not alter the shape but only the ab- 
solute value of the form factor. Thus, the 7c~y transition form factor sets an 
example of a simple exclusive observable for which collinear factorization is 
insufficient for Q 2 as large as 40 GeV 2 . 

It is perhaps of interest to examine the role of the soft intrinsic transverse 
momentum or 6-dependence of the wave function further. Thus, one may 
wonder whether or not a possible scale dependence of the transverse-size 
parameter improves the fit. This scale dependence has been investigated 
in [32] and found that a n slowly decreases with increasing Q 2 . Employing 
the numerical results given in [32] in the fit, one obtains results for the 
form factor that are very similar to the ones presented in Fig. 4. As in most 
applications of light-cone wave functions an approximately scale independent 
transverse-size parameter is therefore assumed in this work for convenience. 
If the transverse-size dependence in (5) is neglected which leads to the original 
version of the MPA [15], a food fit to the form factor data cannot be achieved, 
the results are too flat as compared to the data and the minimal x 2 is 155. 
Hence, the b dependence of the wave function is an important ingredient of 
the MPA as has been suggested in [26]. A 6-dependence of the wave function 
like 

^ P ocexp[-b 2 /4^] (20) 
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leads to fits of similar quality as with (5). The obtained values of the Gegen- 
bauer coefficients are slightly larger and have larger errors. Within these 
errors they however agree with those quoted in (17), (18) within one stan- 
dard deviation. 

One may evaluate the probability (9) of the pion's valence Fock state and 
the r.m.s. value of (8). From the parameters quoted in (17) and (18) one 
finds P q q ~ 0.06 and (k\) l/2 ~ 710 MeV. These values appear to be plausible 
at a low scale although, in view of the assumption of a scale-independent a n , 
they are to be taken with a grain of salt. 

Since the integrations in (1) extend from x — to 1 and from 6 = 
to 1/Aq CD the soft regions contribute to the form factor too. The decisive 
question is how much. In order to investigate this issue a cut-off parameter 
\x c is introduced and any contribution to the form factor is set to zero if the 
renormalization scale (14) is less than /i c . The accumulation profile is then 
defined by the ratio F 7r7 ( / u c )/F 7r7 (0). It is shown in Fig. 5 and reveals that in 
the soft regions defined by small values of the renormalization scale (14), say, 
less than 1 GeV, only small contributions are accumulated. The bulk of the 
contributions to the transition form factor is generated in regions where the 
renormalization scale is sufficiently large. Hence, within the MPA, the 7T7 
transition form factor can be considered as being calculated self-consistently. 
This property holds down to values of the photon virtuality of about 2 GeV 2 . 

Li and Mishima [10] also applied the MPA to the 717 transition form 
factor and achieved a reasonable fit to the CLEO [5] and BaBar data [4]. In 
contrast to the present work the flat distribution amplitude [6] 

= 1 (21) 

is used. It is combined with a Gaussian 6-dependence as in (5) in a kind of 
wave function. However, this product cannot be considered as a proper wave 
function in so far as it is not normalizable, see (9) 3 . It is furthermore argued 
in [10] that the distribution amplitude (21) is accompanied by a threshold 
factor that represents resummed double logs a s In 2 x and a s In 2 (1 — x) arising 
from the end-point singularities which occur for the flat distribution ampli- 
tude in collinear factorization. The threshold factor combined with the flat 

3 The Gegenbauer series of the flat distribution amplitude truncated at uq is however 
normalizable provided no is not too large. Alternatively one may use (20) which leads to 
a normalizable wave function and provides similar results for the transition form factor if 
the transverse size parameter is appropriately chosen. 
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Figure 5: Left: The accumulation profile of the 7T7 transition form factor in 
the space (solid line) and time-like (dashed line) regions versus a cut-off of the 
renormalization scale at Q 2 = 10 GeV 2 evaluated for fit (18). The absolute 
value of the time-like form factor is displayed. (Colors online). Right: The 
effective power of the distribution amplitude (22) compared to the power of 
the threshold factor in [10] at the scale [l. 

distribution amplitude can be viewed as an effective distribution amplitude 
of the type 

•f^W' < 22 > 

According to [10], r is about 1 for low Q 2 and small for Q 2 ~ 40 GeV 2 . This 
particular Q 2 -dependence of the power r generates the increase of the form 
factor required by the BaBar data [4]: At low Q 2 the effective distribution 
amplitude is the asymptotic one implying small values of the transition form 
factor while, at large Q 2 , the effective distribution amplitude is close to the 
flat one and hence leads to much larger values of the form factor. It is to be 
stressed that QCD evolution of the proper distribution amplitude, <& n — 1, 
is omitted in [10]. Since the gliding factorization scale //p is an essential part 
of the MPA this seems to be a problematic assumption. 

Eq. (22) defines a family of power-like distribution amplitudes. It includes 
the limiting cases of the asymptotic distribution amplitude for r = 1 as 
well as the flat distribution amplitude (21) for r = 0. Also the square root 
distribution amplitude proposed in [33] belongs to this family. By comparison 
with the Gegenbauer series (6) of (22) and its evolution one can show [28] 
that the distribution amplitude (22) approximately remains power-like under 
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evolution, r — > r(/i) over a large range of the scale. 

One may examine the power-like distribution amplitude (22) by fitting 
the transverse size parameter as well as the power r(/x ) to the data on the 
7T7 form factor. One finds (/io = 2 GeV) 

a w = 0.40 ±0.05 GeV" 1 r(^ ) = 0.59 ± 0.06 , X 2 = 34.4. (23) 

The quality of this fit to the data on the ttj form factor is similar to that 
presented in [10]. In Fig. 5 the power r(fx) for the fit (23) is compared to the 
scale dependence of the threshold factor used in [10] (in this work the power 
is set to unity for /i 2 <4 GeV 2 ). As can be seen from Fig. 5 the distribu- 
tion amplitude (22) exhibits the usual evolution behavior, it monotonically 
evolves into the asymptotic one, $as, for ji — > oo. On the other hand, the 
scale dependence of the threshold factor or the effective distribution ampli- 
tude advocated for in [10], is drastically different. It is also to be stressed 
that in [10] the threshold factor is evaluated at the scale Q and not at the 
factorization scale /if although it is to be understood as part of the wave 
function. 

The first few Gegenbauer coefficients of the power-like distribution am- 
plitude (23) are compared to those from the fits (17) and (18) in Tab. 1. 
The first Gegenbauer coefficients of the flat and the square root distribution 
amplitude advocated for in [33], are also displayed in Tab. 1 for comparison. 
The comparison is made at a scale of 1 GeV which seems to be a plausible 
value for soft wave functions like (21). While the Gegenbauer coefficients 
of the flat distribution amplitude are generally larger than those quoted in 
(17) or (18), are the coefficients of (23) and the square root distribution 
amplitude smaller. A wave function like \1/ = c exp [— b 2 /4cx 2 ] to which the 
flat distribution amplitude is associated, can be used in a MPA calculation 
without entailing infrared singularities even if evolution is ignored. On the 
other hand, the use of (21) within the collinear factorization approach leads 
to an infrared- singular result which necessitates a regularization prescription 
as for instance the insertion of a mass parameter into the denominator of 
the quark propagator [7]. Such a regularization prescription typically al- 
ters the asymptotic behavior of the scaled form factor drastically: instead of 
Q 2 F 7r7 — > const, the scaled form factor increases logarithmically [6, 7, 8]. 
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04 




(17) 


0.27 ±0.07 








(18) 


0.25 


0.01 ±0.06 





(23) 


0.14 


0.05 


0.03 


(21) 


0.39 


0.24 


0.18 




0.15 


0.06 


0.03 



Table 1: Lowest Gegenbauer coefficients for various distribution amplitude 
at the scale 1 GeV. 

4 Generalization to the rj and 7/7 form factors 

The analysis of the 7T7 transition form factor within the MPA can straight- 
forwardedly be extended to the cases of the 777 and 7/7 ones [34, 35]. These 
transition form factors may be expressed as a sum of the flavor-octet and 
flavor-singlet contributions (P = 77, r/) 

F Pl = F% + Fp 7 . (24) 

As is the case for the form factor the functions F Pl (i = 1,8) are pro- 
portional to the constants f P assigned to the decays of meson P through the 
SU(3)p octet or singlet axial- vector weak currents which are defined by the 
matrix elements 

(0 I JU I P(P)) = ifpP, (25) 
Adopting the general parameterization [36] 

/« = / 8 cos^ 8 , fl = -/xsinflx, 

f} = issuing, $ = A cos 0i, (26) 

one can show [37] that on exploiting the divergences of the axial-vector cur- 
rents - which embody the axial-vector anomaly - the mixing angles, 8s and 
61, differ considerably from each other and from the state mixing angle, 8. 
In [37] the mixing parameters have been determined: 

/ 8 = 1.26 U, h = 1-17 U, 

9 8 = -21.2°, e l = -9.2°. (27) 
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Assuming particle-independent [34, 37] octet and singlet wave functions, 
\l/ 8 and ty 1 , which are parameterized as in (5) with the respective decay 
constants f$ and f\ instead of f n , one can cast the transition form factors 
into the form 

F m = cos # 8 -F 8 — sin 6i F 1 

F^ 7 = sin # 8 F 8 + cos 9i F 1 . (28) 
The charge factors in (3) read (with P = 1,8) 

Cs = (e 2 u + e 2 d - 2e 2 s )/V6, C x = (e 2 u + e 2 + e 2 s )/V3. (29) 
The asymptotic behavior of the form factors is 

Q^^f-fs, Q 2 F^^=h. (30) 

It is to be noted that the singlet-decay constant is renormalization-scale 
dependent [36]: 

^fu = 7A(/i) fi (31) 
where the anomalous dimension is 

1A = -n f (a s (fi)/nf . (32) 

Since the anomalous dimension controlling this scale dependence is of order 
a 2 , it leads to tiny effects. In fact, if the value of fi quoted in (27) is to be 
understood as being valid at, say, the scale ~ 1 GeV, its asymptotic value is 
13% smaller; over the range of available data it only decreases by less than 
2%. The scale dependence of f\ is therefore discarded for convenience. 

It is to be stressed that to NLO of the hard scattering there is also a 
contribution from the glue-glue Fock component of the mesons to the singlet 
form factor [38, 39]. In the MPA analysis performed in this work the glue-glue 
Fock component does not contribute directly but only through the matrix of 
the anomalous dimensions and the mixing of the singlet quark- ant iquark with 
the glue-glue distribution amplitude. It is assumed here that the Gegenbauer 
coefficients of the glue-glue distribution amplitude are zero at a low scale 
of order 1 GeV. Hence, the quark- ant iquark distribution amplitude evolves 
with the eigenvalues 7^ + ) of the anomalous dimension matrix [39] . The values 
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Figure 6: The octet and singlet form factors. Dotted lines represent the 
asymptotic behavior (30), the dashed lines the results obtained in [34]. The 
solid lines represent the new fit (33). Data taken from [5, 19]. (Colors online) 

of 7^ + ) practically fall together with those of 7„, the anomalous dimension 
of the octet distribution amplitude, see (7). For instance, = 0.639 while 
72 = 2/3 for four flavors. 

The two form factors F 8 and F 1 can now be evaluated from (1) and (3) 
in full analogy to the transition form factor. The data on F 8 and F l 
are extracted from the CLEO [5] and BaBar [19] data using (28). As for 
the 7T7 form factor only the transverse size parameter and one Gegenbauer 
coefficient for each wave function can be determined. The best fit is obtained 
with the parameters (fio = 2 GeV): 

a 8 = 0.84 ±0.14 GeV" 1 , a 8 2 (fx ) = -0.06 ± 0.06 , 

o-i = 0.74 ±0.05 GeV" 1 , a\(n Q ) = -0.07 ±0.04. (33) 

The values of \ 2 are 15.0 and 14.1 for the octet and singlet cases, respectively 
(for 16 data points in each case). The probabilities of the singlet and octet 
wave functions are 0.24 and 0.19, respectivly. The corresponding r.m.s. k^- 
values are 390 and 440 MeV. 

In Fig. 6 the results of this fit are compared to the data on F 8 and F 1 . 
The quality of this fit is very good. In contrast to the ttj case the data on 
both F 8 and F 1 lie below the asymptotic results (30). The combination of 
these two form factors into the physical ones according to (28) leads to the 
results which are shown in Fig. 7. Again very good agreement with the data 
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Figure 7: The scaled 777 and 7/7 transition form factor versus Q 2 . Data taken 
from [5, 18, 19]. For notations refer to Fig. 6. (Colors online) 

is to be observed. Also shown in Figs. 6 and 7 are the results obtained in 
[34] which have been evaluated from the asymptotic distribution amplitudes 
(with <Ti = cr 8 = 0.861 GeV -1 ). The octet as well as the 7/7 form factors of 
[34] are in very good agreement with experiment while the results for F 1 and 
the 7/7 form factor are somewhat too large. 

The form factors scaled by their respective asymptotic behaviors are dis- 
played in Fig. 8. All the form factor ratios for the light mesons behave 
similar although not identical. Asymptotically they all tend to 1. The irj 
form factor approaches 1 from above, the other ones from below. The ap- 
proach to 1 is a very slow process; even at 500 GeV 2 the limiting behavior 
has not yet been reached. It is also evident from Fig. 8 that, forced by 
the BaBar data, there are strong violations of SU(3)i? flavor symmetry in 
the groundstate octet of the pseudoscalar mesons at large Q 2 . The differ- 
ence between the n / ~f and the 777 or more precisely the 77s7 form factors is 
larger than the difference between their respective decay constants. In other 
processes involving pseudoscalar mesons, e.g. two-photon annihilations into 
pairs of pseudoscalar mesons [40, 41], such large flavor symmetry violations 
have not been observed. Below 8 GeV 2 , i.e. in the range of the CLEO data, 
flavor symmetry breaking is much milder. The r/ c 7 transition form factor 
which is also shown in Fig. 8, behaves very different. This form factor will 
be discussed in the next section. 
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Figure 8: Left: The P — 7 transition form factors scaled by the corresponding 
asymptotic behavior, versus Q 2 . The thick solid (dashed, dotted, thin solid) 
line represents the case of the 77 (r/, ir, t] c ) meson. Parameters are taken 
from fit (18), (33) and (37). Right: The 7/ c 7 form factor scaled by its value 
at Q 2 = 0. Data taken from [20]. The solid (dotted) line represents the 
results of a calculation with the values of the parameters quoted in (37) 
(with m c = 1.26 GeV). The dashed line is the prediction given in [42]. 
(Colors online) 



5 The case of the r\ c 

The essential difference of the 7^7 transition form factor to the other three 
form factors is the large mass of the 7] c (M Vc ) or that of the charm quark (m c ). 
They constitute a second large scale in addition to the virtuality of one of 
the photons. The masses cannot be neglected in a perturbative calculation in 
constrast to the case of the light mesons where quark and hadron masses do 
not play a role. Despite this the r/ c 7 form factor is to be calculated from (1) 
but the hard scattering amplitude to lowest order perturbative QCD reads 
[42] 



T = ^ (M) 

H xQ 2 + (1 + Axx)m 2 c + k 2 ± ' 1 ; 

The symmetry of the problem under the replacement of x by x is already 
taken into account in (34). Due to the involved second large scale in the 
problem the ^7 form factor can be calculated even at Q 2 = 0. 
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The light-cone wave function of the r] c is parameterized as in (5) 4 . Fol- 
lowing [42, 43] the distribution amplitude is chosen as 



$r, c = N(cr Vc )xxexp - a 2 M, 



2 (x - 1/2) 2 



(35) 



XX 



where N(a Vc ) is determined from the usual requirement J dx$ Vc (x) = 1. 
The distribution amplitude exhibits a pronounced maximum at x — 1/2 and 
is exponentially damped in the end-point regions. It describes an essentially 
non-relativistic cc bound state; quark and antiquark approximately share 
the meson's momentum equally. In the hard scattering amplitude the charm 
quark mass occurs while in the distribution amplitude the meson mass is 
used. This property of the latter distribution amplitude is a model assump- 
tion which contributes to the theoretical uncertainty of the results. In the 
sense of the non-relativistic QCD [44] 2m c and M Vc are equivalent. In (1) 
the Sudakov factor exp [-S] can be set to 1 in the case at hand for two rea- 
sons: First, due to the large, non-negligible c-quark mass the radiative QCD 
corrections only produce soft divergencies but no collinear ones and, hence, 
the double logs do not appear. Second, the Sudakov suppressions is mainly 
active in the end-point regions (c.f. the discussion in Sect. 2) which are al- 
ready strongly damped by the f] c wave function. The evolution behavior of 
the 7] c distribution amplitude is unknown in the range where Q 2 is of order 
of M 2 c and is therefore ignored here. Consequently, also the running of the 
charm quark mass is omitted. It has been checked that the effect of the this 
scale dependence is anyway only on the percent level. 

The normalization of the 7^7 transition form factor is fixed by its value 
at Q 2 = which is related to two-photon decay width by 



However, this decay width is experimentally not well known [48] . It is there- 
fore advisable to normalize the form factor by its value at Q 2 = all the 
more so since the recent BaBar data [20] are also presented this way. Do- 
ing so the perturbative QCD corrections at Q 2 = to the r/^ transition 
form factor which are known to be large [45], are automatically included. 

4 In [42] the Gaussian (20) is taken. The version (5) is chosen here in order to be 
conform with the calculations of the other form factors. The differences between the two 
versions are marginal as has been mentioned previously. 



r[^ c ->■ 77] 



-na 2 lm Ml 



^ c7 (0) 



(36) 
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Also the a s corrections for Q 2 < M 2 [46] cancel to a high degree in the ratio 
F r)c7 (g 2 )/F r)c7 (0). Even at Q 2 = 10 GeV 2 their effect is less than 5%, c.f. 
the discussion in [42]. The uncertainties in the present knowledge of the r] c 
decay constant do also not enter the predictions for this ratio. 

The recent Babar data on the r/^ form factor [20] are shown in the right 
hand panel of Fig. 8. The behavior of this data has indeed been predicted 
in [42]. The predictions which have been evaluated from m c = M Vc /2, are 
about one standard deviation too large but with regard to the uncertainties 
of the theoretical calculation, as for instance the exact value of the mass of 
the charm quark, one can claim reasonable agreement between theory and 
experiment. A little readjustment of the value of the charm quark mass 
improves the fit. Thus, with the parameters 

m c = 1.35 GeV, a Vc = 0.44 GeV" 1 , (37) 

a perfect agreement with experiment is achieved as is to be seen in Fig. 8. For 
comparison there are also shown results evaluated from m c = 1.21 GeV in 
Fig. 8. As one may note from the left hand panel of Fig. 8 the 7^7 transition 
form factor behaves quite differently from the other three form factors. The 
large charm-quark mass slows down the approach to the asymptotic limit 

Q 2 F Vcl -> ^ • (38) 

Finally one may examine whether the parameters quoted in (37) are plau- 
sible. For this purpose the r/ c 7 form factor at zero momentum transfer is 
evaluated. Ignoring a s -corrections as well as relativistic effects and using the 
value 420 MeV for the r/ c decay constant f Vc , one obtains 

F„ c7 (0) = 0.085 GeV" 1 , (39) 

and a two-photon decay width (36) of 

r(?7 c ->■ 77) = 8.05 keV . (40) 

This result is in good agreement with the value of (7.20±2.11) keV evaluated 
by the PDG [48] and with recent theoretical estimates, see for instance [47] 
and references therein. The parameters quoted in (37) together with f Vc = 
420 MeV correspond to a normalization of the r\ c distribution amplitude 
N(a Vc ) = 8.849, to a probability of the valence Fock state of 0.82 and to 
a r.m.s. of 773 MeV. The latter two values appear reasonable for a 
quarkonium state. 
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6 Two virtual photons 



An extension of the MPA analysis of the P7 transition form factors to the 
case of two virtual photons can be found in [49, 50, 51]. As an example the 
7T7* form factor will be discussed here. Denoting the photon virtualities by 
Q 2 and Q' 2 and introducing the variables 

<?^« 2 + Q' 2 ), «=$^. («) 

one finds for the hard scattering amplitude the expression 

f H = -?L-K {yJl-u(l-2x)Qb) (42) 

V37T 

in 5-space. In generalization of (1) the nj* form factor now reads 
F nr (Q 2 ,uj) = y^^^(z,-b,// F )T„(z,^ 

(43) 

The renormalization scale is taken as 

H R = max (y/l - u(l - 2x)Q, 1/6) (44) 
in generalization of (14). 

2 2 

The form factor only falls off like 1/Q at large Q in contrast to the 
Q~ 2 Q'~ 2 oc Q~ 4 behavior of the vector meson dominance model [52]. It 
is also interesting to note that like in collinear factorization, the 717* form 
factor is insensitive to the higher Gegenbauer terms for Q 2 ~ Q' 2 ; in the limit 
uj — > it only depends on the asymptotic distribution amplitude. In fact it 
has been shown [50, 51] that the n-th order Gegenbauer term contributes 
to order u n . The importance of effects from the intrinsic kj_ diminishes as 
both photons become virtual. This is evident from the quark propagator in 
momentum space reading 

T H oc -2x))Q 2 + k 2 ± }- 1 . (45) 

As u deviates from 1, the real photon limit, the form factor F ny * becomes 
less sensitive to the end-point regions where either the quark or the antiquark 
becomes soft. In the limit 00 — > (45) reduces to [Q + k^] _1 . This is to 
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be contrasted with the ou — > 1 limit in which T H oc [xQ 2 + k 2 ] 1 , see (2). 
Thus, kj_ plays a minor role in the quark propagator provided u is sufficiently 

— 2 

small. In the limit u — > and large Q the form factor becomes 

Q 2 F nr (Q 2 ,u) = ^f n , (46) 

a result that has been derived long ago in [53] . 

The measurement of the space-like 717 transition form factor is performed 
in e + e~ collisions with so-called single-tag events where either the electron 
or the positron in the final state is detected. This method implies an integra- 
tion over the spectrum of that exchanged photon which is emitted from the 
undetected lepton, up to a value of, say, Q' 2 set by an experimental cut. Ac- 
tually for the BaBar experiment the cut is Q' 2 < 0.18 GeV 2 . Thus, strictly 
speaking, one measures the form factor for the transition from a quasi-real 
photon with an effective virtuality Q' e g less than the cut value, to the pion. 
The question arises how good does this measured form factor approximate 
the one for the transition from a real photon to the pion. In order to examine 
this issue F nj * is evaluated for small values of the ratio Q' 2 /Q 2 using (43) and 
the pion wave function parameters from fit (18). As one may see from Fig. 
9 where the results for a sample set of Q 2 values are shown, the transition 
form factor depends on Q' 2 mildly (at least within the MPA) as long as the 
ratio Q' 2 /Q 2 is smaller than about 0.01. At this value of the ratio of the two 
virtualities the form factor is reduced by about 5% as compared to its value 
at the real photon limit while at a value of 0.1 of the ratio, the form factor 
is smaller by about 30%. 

7 Remarks on the time-like transition form 
factors 

In the context of the large Q 2 behavior of the transition form factors one 
may also consider the recent measurement of the time-like 777 and 7/7 form 
factors at s = 112 GeV 2 by the BaBar collaboration [54]: 

s\F V7 \ = 0.229 ±0.031 GeV, s|F^ 7 | = 0.251 ± 0.021 GeV. (47) 

Comparing these data with the predictions for the space-like form factors 
which are shown in Figs. 7 and 8, one notices that for the case of the rj the 
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Figure 9: The ir^y* form factor versus the ratio of the two photon virtualities 
for Q 2 = 25(10,4) GeV 2 represented as solid (dashed, dotted) line. The 
parameters of the pion wave function are taken from fit (18). (Colors online) 

experimental time-like value is much larger than the predicted space-like one 
which amounts to 0.169 at Q 2 = 112 GeV 2 . For the case of the r)', on the 
other hand, the experimental time-like value practically falls together with 
theoretical result of 0.258 for the space-like form factor. 

A detailed study of the pseudoscalar meson-photon transition form fac- 
tors in the time-like region is beyond the scope of the present article. Within 
the MPA however, one expects slightly different values for the time-like and 
space-like transition form factors: The space- like propagator [xQ 2 + k^] -1 
is to be replaced by [— xs + — it]" 1 in the time-like region or, in 6-space, 
K (^Qb) in (3) by i7r/2H^\y/xsb) with being the zeroth order Han- 
kel function. The pole of the propagator now occurs within the range of 
integration and in general leads to an enhancement as well as a phase of the 
form factor. For the case of the electromagnetic form factor of the pion this 
phenomenon has been pointed out and studied in some detail by Gousset and 
Pire [55]. The analytic continuation of the Sudakov factor from the space- 
like to the time-like region, which is necessary too, is not well understood. 
It probably leads to an oscillating phase [55, 56]. 

In order to give an admittedly rough estimate of the expected size of 
the time-like form factor it is followed [55] and the time-like propagator is 
combined with the space-like Sudakov factor. With this receipe one finds 
the absolute values 0.25 and 0.24 GeV for the scaled time-like 717 transition 
form factor at s = 3 and 100 GeV 2 , respectively. The corresponding ratios 
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of the time and space- like transition form factors are ~ 1.9 and ~ 1.1. For 
s larger than about 5 GeV 2 the time-like form factor is dominantly real, 
its imaginary part contributes less than about 10% to the absolute value. 
The accumulation profile of the absolute value of the time-like form factor 
is displayed in Fig. 5. Although the pole of the propagator occurs in the 
end-point region (at x = k±/s) the profile is only mildly softer than that 
of the space-like form factor, i.e. also the time-like form factor is mainly fed 
by contributions from regions where the renormalization scale is sufficiently 
large. In the light of this feature one may consider the results for the time-like 
transition form factor obtained within the MPA as a tolerable estimate. 

Along these lines one can also compute the 777- and 7/7 transition form 
factors in the time-like region. Using the parameters quoted in (33) one 
finds that at s = 3 GeV 2 the absolute values of the time-like form factors are 
about a factor of 2 larger than the space-like form factors. At s = 112 GeV 2 
the results are s|F^ 7 | ~ 0.17 GeV and s|F^ 7 | ~ 0.28 GeV. While for the 
case of the 7/ there is rough agreement with experiment (47) within errors, 
is the 777 form factor too small by about two standard deviations. 

8 Concluding remarks 

In this paper an analysis of the form factors for the transitions from a photon 
to a pseudoscalar meson is presented. The analysis is performed within the 
MPA which bases on factorization. It is shown that due to the Sudakov 
suppressions which are an important ingredient of the MPA and which repre- 
sents radiative corrections in next-to-leading-log approximation summed to 
all orders of perturbation theory, higher order Gegenbauer terms of the me- 
son's distribution amplitude are suppressed at low Q 2 . In fact, the combined 
effect of the hard scattering kernel and the Sudakov factor leads to a series of 
power suppressed terms which stem from the soft regions. The intrinsic k^- 
dependence generates a second series of power corrections which in contrast 
to the first series, do not grow with the Gegenbauer index. The interplay 
of these two ingredients results in a remarkable feature of the MPA - the 
transition form factors are only affected by the few lowest Gegenbauer terms 
of the distribution amplitude, the higher ones do practically not contribute. 
How many Gegenbauer terms are relevant depends on the range of Q 2 consid- 
ered: In the Q 2 range covered by the CLEO data [5] (< 10 GeV 2 ) it suffices 
to use just the asymptotic distribution amplitude in order to fit the CLEO 



26 



data. With the BaBar data [4, 19] at disposal, covering the unprecedented 
large range 4 GeV 2 < Q 2 < 35 GeV 2 , the next or the next two Gegenbauer 
terms have to be taken into account or, turning the argument around, can 
be determined from an analysis of the data on the transition form factors. 
Indeed this is what has been done in this work. From the present analysis 
it turns out that for the case of the pion a fairly strong contribution from 
a2 is required by the data while for the i] and rj' much smaller deviations 
from the asymptotic distribution amplitude are needed. For these cases the 
results from a previous calculation within the MPA [34] are already in fair 
agreement with the BaBar data, nearly perfect for the rj, slightly worse for 
the rf . Comparing the 717 form factor with the 777 or more precisely the 
7^7 one, one observes a strong breaking of flavor symmetry in the ground- 
state octet of the pseudoscalar mesons. The difference between the two form 
factors is larger than the respective decay constants. In other processes in- 
volving pseudoscalar mesons such large flavor symmetry violations have not 
been observed. With regard to the theoretical importance of the transition 
form factors, in particular the role of collinear factorization a remeasurment, 
e.g. by the BELLE collaboration, would be highly welcome. 

One may wonder what the implications of the new distribution amplitudes 
for the pseudoscalar mesons are for other hard exclusive processes. A detailed 
investigation of this issue is beyond the scope of the present paper. It has 
however been checked that for the pion's electromagentic form factor there 
is no substantial change. The perturbative contribution to it still amounts 
to only about a third of the experimental value of the form factor which 
is measured only at low Q 2 [62]. The perturbative contribution is slightly 
increasing with growing Q 2 now and not as flat as in [26] . In any case for Q 2 
less than 10 GeV 2 the differences are marginal. 
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A Appendix 



In this appendix details of the Sudakov factor are presented. The Sudakov 
exponent S which comprises the characteristic double logarithms produced 
by overlapping collinear and soft divergencies (for massless quarks) has been 
calculated in [14]. In axial gauge, n ■ A — 0, these overlaps arise in general 
from two-particle reducible Feynman graphs where, before the hard interac- 
tion, the gluon is exchanged between the quark and antiquark of the meson 
or emitted from and reabsorbed by either the quark or the antiquark. The 
Sudakov factor can therefore be analyzed indepedently of the physical pro- 
cess and can be viewed as part of the meson wave function. For the case of 
interest, only the two-particle reducible graphs occur anyway. As shown in 
[14] for a quark-antiquark system the Sudakov exponent reads 5 

S = s(x, b, Q) + s(l -x,b,Q) + 2 [ m ^ >y q (a B (fi)) • (A.l) 

The impact parameter b, canonically conjugated to k^, is the spacial sep- 
aration of quark and antiquark. The Sudakov function s(£, b, Q) where £ is 
either x or 1 — x, is given by 

s(Z,b,Q) = \ 1°** — \2ln(^-)A(C 1 ,g( f i)) + B(C 1 ,C 2 ,g^))\ . 

(A - 2) 

An appropriate choice of the gauge vector n has been made in order to obtain 
this result. The function B has been calculated to order a s in the MS scheme 
explicitly. The function A arises from the use of the renormalization group 
equation which is applied in order to absorb all the Q 2 dependence into the 
scale of the coupling constant. It has been calculated to order a^: 

A= ^ A (i) + (^L) 2 A m + (D (a 3 s ) (A.3) 

7T V 7T / 



with 



2 



67 7T 2 ^ 

-- Y --n f + f3oln(C ie ^/2) 



(A.4) 



5 Musatov and Radyushkin [57] claim that, due to kinematical properties the Sudakov 
functions should be summed as s(x, 6, Q) + s(\/E, b, Q) where it is assumed that the real 
photon is attached to the antiquark line. Doing so the universality of the wave function 
(including the Sudakov factor) is broken. This alternative possibility only leads to tiny 
numerical differences in the form factor. The predictions change by less than 1%. 
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where 7^ is the Euler constant. The constants C\ and C2 are free parameters 
which may be chosen in such a way that large logarithms from higher orders 
in the perturbative expansion of B are avoided. The conventional choice is 

y/2. 



Ci 



Co 



(A.5) 



Other choices of the Cj have been discussed in [58]. 

The Sudakov function has been explicitly given in [14] first. Later rep- 
etitions of this calculation [59, 60] led to a slightly different result which is 
quoted here 



s(£,b,Q) = 



+ 



+ 4: 
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PI 
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ln(2a + l_ M 2S) + l ) + l ln2(2 , ) _l ln2(2S) 
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The variables q and 6 are defined by 



lnf- 
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(A.6) 



(A.7) 



"2A QCD y v 6Aq CD - 

The last term in (A.6) represents the integrated function B. The differences 
between the various results for s used earlier [14, 15, 26] and (A.6) are nu- 
merically small, of the order of a few percent. It should be noted that the 
last term in (A.6) is twice as large as in [58]. Also this discrepancy has little 
effect on the form factor. 

The integral in (A.l) arises from the application of the renormalization 
group equation. It combines the effects of the renormalization group equa- 
tion on the wave function (with the factorization scale (ip) and on the hard 
scattering amplitude involving the renormalization scale /ir. The evolution 
from one scale to another is controlled by the anomalous dimension of the 
quark's wave function in axial gauge [61] 7 g = —a s /n + C(o; 2 ). It leads to 



J Hp 



7«(«b(a*)) 



2_ ^ lnj/4/A|cD) 



(A.8) 



» ' Po"' In (/4/ A qcd) 

with the help of the 1-loop result for a s . Its use is consistent with the deriva- 
tion of (A.6). In particular the terms in the first line of (A.6) are calculated 
from the 1-loop a s and these terms dominate the Sudakov function. 
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